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Abstract. Loop quantum cosmology applies techniques derived for a background independent 
quantization of general relativity to cosmological situations and draws conclusions for the very 
early universe. Direct implications for the singularity problem as well as phenomenology in the 
context of inflation or bouncing universes result, which will be reviewed here. The discussion 
focuses on recent new results for structure formation and generalizations of the methods. 


1. Introduction 

The distinguishing feature of general relativity, in comparison to other interactions, is the fact 
that the metric as its basic field does not just provide a stage for other fields but is dynamical 
itself. In particular in cosmological situations the metric differs significantly from a static 
background and cannot be written as a perturbation. Thus, a faithful quantization requires 
a background independent formalism which then must be non-perturbative. 

An approach to quantum gravity which realizes this from the outset is loop quantum gravity 
dnnsm. Here, background independence leads to a discrete structure of geometry whose 
scale is a priori free (set by the Barbero-Immirzi parameter 0 El; in this paper we set the 
value equal to one for simplicity) but fixed to be close to the Planck scale by black hole entropy 
calculations 0 Ei ei uni- Thus, it is not of relevance on directly accessible scales and will only 
become noticeable in high curvature regimes. In particular, this is the case close to the big bang 
where the universe itself is small. 

Classically, the universe would emerge from or evolve into a singularity at those scales, where 
energy densities blow up and Einstein’s equations break down. For a long time, it has been 
hoped that quantum gravity will resolve this problem and provide a more complete framework 
which does not break down. Moreover, since this will inevitably come with modifications of 
the classical theory at small scales, one can expect phenomenological and potentially observable 
consequences in the very early universe. 

Even classically, it is difficult to analyze the situation in full generality, and the quantum 
theory is even more complicated and less understood. A common strategy in such a situation 
consists in introducing symmetries which can be taken as homogeneity or isotropy in the 
cosmological context. In contrast to earlier approaches initiated by Wheeler and DeWitt mm, 
the theory has now been developed to such a level that the introduction of symmetries can be 
done at the quantum level by employing symmetric states m, rather than reducing the classical 
theory first and then quantizing. The relation to the full theory is thus known, and it is possible 
to ensure that special features required for a consistent background independent formulation 
translate to the symmetric context. 


It is then possible to take properties of the full theory, transfer them to symmetric models 
and analyze them in this simpler context. In particular, the discreteness of spatial geometry 
survives the reduction M, which is already a difference to the Wheeler-DeWitt quantization. 
It also implies that there are in fact modifications at small scales coming from the full theory, 
whose phenomenological consequences can be studied in cosmological models m- 


2. Variables 

A spatially isotropic space-time has the metric 

ds 2 = -df 2 + (1 ^ r2)2 dr 2 + a(f)Vdfl 2 

where k can be zero or ±1 and specifies the intrinsic curvature of space, while the scale factor 
a(t) describes the expansion or contraction of space in time. It is subject to the Friedmann 
equation 

3(a T k)a — 8TrGH mSL n eT {^(i, (1) 

where G is the gravitational constant and iL ma tter the matter Hamiltonian (assumed here to be 
given only by a scalar cj) and its momentum p^). The matter Hamiltonian depends on the matter 
fields, but also on the scale factor since matter couples to geometry. In the case of a scalar, for 
instance, we have 

-^matter = + a 3 ^ {(j>) (2) 

with the scalar potential V(cf>) and the classical momentum = a 3 (j). 

Loop quantum gravity is based on Ashtekar variables, which provide a canonical formulation 
of general relativity in terms of a densitized triad and an SU(2) connection on space. In 
the isotropic context this reduces to working with the isotropic triad component p with 
\p\ = a? and the isotropic connection component c = ^(k + d) which are canonically conjugate: 
{c,p} = 8irG/3. There is one essential difference to the metric formulation: p can take both 
signs since it depends on the orientation of the triad. Thus, p does not only determine the 
size of space through \p\, but also its orientation via sgnp. (Another difference is that, when 
isotropic models are derived through homogeneous ones, a canonical formalism is not available 
for k = — 1. We will thus restrict ourselves to k = 0 and k = 1.) 

Dynamics in the canonical formulation is dictated by the Hamiltonian constraint 

H = —3(47rG)' 1 [2c(c - k) + k 2 ] y/\p\ + H matteT (p, ^>,^) = 0 (3) 


which indeed reduces to the Friedmann equation upon using the definition of p and c. Moreover, 
the Hamiltonian constraint H gives Hamiltonian equations of motion for gravitational variables, 
such as c = {c, H} resulting in the Raychaudhuri equation 
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( 4 ) 


and matter equations of motion, e.g. for a scalar 

<t> = {4>,H} = p^/a 3 
V4> = {P<t» H } = -o 3 R'(</>) 


which lead to the Klein-Gordon equation 

4> + 3aa~ l (j) + V'{(j)) = 0 . 


( 5 ) 





3. Loop quantization 

While a Wheeler-DeWitt quantization would start with a Schrodinger representation and work 
with wave functions ip(a) such that a is represented as a multiplication operator and its 
momentum related to a by a derivative operator, the loop quantization implies an inequivalent 
representation m- Here, one usually starts in the connection representation such that states 
are functions of c, an orthonormal basis of which is given by 

(c\p) = e iflc/2 , fjtsR. (6) 

Since these states are by definition normalized, it is clear that the Hilbert space is non-separable 
(it does not have a countable basis) and that the representation is inequivalent to that assumed 
in the Wheeler-DeWitt quantization. 

Basic operators, which quantize p and c, also have properties different from a as a 
multiplication operator or its conjugate as a derivative operator. The action of basic operators 
on states © is given by 


p\v) = 

( 7 ) 

e ^' c / 2 \p) = \p + p!) 

( 8 ) 


with the Planck length = s/SuGh. Thus, since all eigenstates \p) of p are normalizable, p 
has a discrete spectrum. Moreover, there is only an operator for the exponential of c, not c 
directly. Both properties are very different from the corresponding operators in the Wheeler- 
DeWitt quantization where the analog of p , the scale factor a, has a continuous spectrum and 
its momentum has a direct quantization. 

On the other hand, the properties of the basic operators 0, © are analogous to those 
in the full theory, where also flux operators quantizing the triad have discrete spectra and 
only holonomies of the connection are well-defined operators but not the connection itself. In 
the full theory, these properties are consequences of the background independent formulation: 
One has to smear the basic fields given by the connection and the densitized triad in order 
to have a well-defined Poisson algebra to represent on a Hilbert space. In field theory this is 
usually done in a three-dimensional way using the background metric to provide a measure. 
This is certainly impossible in a background independent formulation, but there are natural, 
background independent smearings of the connection along one-dimensional curves and of the 
densitized triad along surfaces. Their algebra, the holonomy-flux algebra, is well-defined and 
one can then look for representations. Here, it turns out that there is a unique one carrying a 
unitary action of the diffeomorphism group In this representation, fluxes and 

thus spatial geometric operators which are built from triad components have discrete spectra as 
a direct consequence of background independence. 

This representation is then carried over to symmetric models such that triads have discrete 
spectra, too, and only exponentials of connection components are directly represented. These 
properties of the loop representation define the structure of the algebra of basic operators, and 
they have far-reaching consequences: 
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As a consequence of the discrete triad spectrum, operators quantizing the inverse volume 
are finite despite the classical divergence. This already signals a more regular behavior at the 
classical singularity which has to be confirmed by using the quantum dynamics. This dynamics, 
as a consequence of the second basic quantity that only exponentials of c are represented, 
happens in discrete internal time. Together with the properties of inverse volume operators 
this combines to a non-singular cosmological evolution. Moreover, inverse volume operators 
imply non-perturbative modifications to the classical Friedmann equation, while the second 
basic property leads to perturbative higher order terms. Both corrections have phenomenological 
consequences. 


4. Non-singular evolution 

In this section we discuss the consequences of basic loop properties as for their implication of the 
quantum evolution [22123 EH- In the following section we will then turn to phenomenological 
consequences ESI- 

4-1- Finite inverse volume 

Since, as one of the basic loop effects, the triad operator p has a discrete spectrum containing 
zero it does not have a densely defined inverse. On the other hand, if we want to quantize a 
matter Hamiltonian such as © which enters the dynamics, we always need inverse powers of the 
scale factor in the kinetic term. It seems that quantum cosmology based on a loop representation 
would already come to an end at this basic step. However, there are general methods in the full 
theory m which allow us to find a well-defined quantization of a 3 . To that end we first rewrite 
the classical expression in an equivalent way which is better suited to quantization. Since such a 
rewriting can be done in many classically equivalent ways, this in general leads to quantization 
ambiguities in non-basic operators. For instance, the inverse volume can be written as 

/ 3 \ 3 /( 2 - a ) 

m := a ~ 3 = \ i*GW + m + 1) j (9) 

where j E (denoting the SU(2) representation in which we take the trace of holonomies 
hi = exp (ctj) with SU(2) generators 77 = —iaj/2 in terms of Pauli matrices 07 ) and 0 < l < 1 
are ambiguity parameters. The advantage of these new expressions is that we now have only 
positive powers of p on the right hand side which, as well as the holonomies, we can easily 
quantize. The Poisson bracket will then be turned into a commutator at the quantum level 
resulting in a well-defined operator whose eigenvalues 

\ 3/ (2—21) 

k\Pn+2k\ l J (10) 

on eigenstates \p) follow from the action of basic operators. 

Since this operator is finite E3, the classical divergence of a 3 is now indeed removed as 
can be seen from the eigenvalues. In particular for larger j , which are sometimes helpful for 
phenomenological purposes, the exact expression can be difficult to use and the approximation 
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is helpful. Even though there are quantization ambiguities, the important properties such as 
the finiteness of the operator and its classical limit are robust. 

4-2. Difference equation 

In order to consider dynamics and to decide whether or not the classical singularity persists as 
a boundary to the evolution, we need to quantize the Friedmann equation m- This is most 
conveniently expressed in the triad representation given by coefficients in an expansion 

\ip) = ^ in triad eigenstates. Since we have to use the basic operator © which is 

a shift operator on triad eigenstates, the quantized Friedmann equation becomes a difference 
equation for 

(^+5 - ^+ 3 )^+ 4 (</>) - (2 + k 2 )(V^ +1 - (13) 

+ (Vu-3 - Vfj.- 5 )e~ lk f> fl -4((j)) = -^7rG^p^matter(M)^/i(0) 

in terms of volume eigenvalues = (£p|/r|/ 6 ) 3//2 . There are also possible ambiguities in this 
constraint, for instance analogous to the parameter j above which have been analyzed in ED- 
Moreover, a symmetrized version is possible, which we do not discuss here for simplicity. 

The evolution dictated by this difference equation in internal time /r does not stop at any 
finite value of /r. In particular, we can uniquely evolve initial values for the wave function 
through the classical singularity situated at [i = 0. Thus, there is no singularity where energy 
densities would diverge or the evolution would stop. This comes about as a consequence of 
the basic loop properties: the discreteness of spatial geometry leads to finite operators for the 
inverse volume as well as evolution in discrete internal time. Both properties enter in the 
demonstration of singularity free evolution. Physically, this means that around the classical 
singularity continuous space-time and with it the classical theory dissolve. Discrete quantum 
geometry, on the other hand, still makes sense and allows us to evolve to the other side of the 
classical singularity. 

5. Phenomenology 

The density © does not just give us a kinematical hint for the removal of classical singularities, 
it is also important as an ingredient in matter Hamiltonians such as 0 . Since at small scales 
the classical a~ 3 is modified by CD), we obtain modified Hamiltonian equations of motion and 
a modified Friedmann equation. For a scalar they are the effective Friedmann equation 

3ad 2 = 8nG {\d(a) eS p 2 (j) + a 3 V(4>)) , (14) 

and Raychaudhuri equation 

\ = ~ (o-U(o)iV 2 (l - l tt dlog (^ 1 “W ) ) - v M ) (15) 

for the scale factor and the effective Klein-Gordon equation 

7 } ■ dlogd(a) e ff 3 ,/ \ t r/ / i\ nr\ 

9 = (pa - - - a a(a) e gV {(p) (16) 

da 

for the scalar. Other matter types have been discussed in ED- These modifications from d(a) at 
small scales lead to diverse effects which give us a new picture of the very early universe. Before 
discussing these effects we note that even though small scales behave very differently from large 
ones, there is an interesting duality in the effective equations which can be helpful in analyzing 
solutions m • 







Figure 1. The effective density (left) implying early inflation (right), which is not realized with 
the classical density (dotted). 


5.1. Inflation 

At small a, the effective density d(a) e g ~ a 3 /( 1- 0 is increasing since 0 < l < 1. Thus, in contrast 
to the classically decreasing a -3 the effective density implies a matter energy on the right hand 
side of the Friedmann equation m which increases with the scale factor. Since the negative 
change of energy with volume defines pressure, this quantum geometry effect naturally implies 
an inflationary phase in early stages IS,- As demonstrated in Fig. |T| inflation automatically 
ends when the peak of the effective density is reached such that there is no graceful exit problem. 

Since the modification is present in the kinetic term of any matter field, we do not need to 
assume <j) to be an inflaton with special properties. Thus, there are several different scenarios 
depending on whether or not we assume an inflaton field to drive the expansion. It is, of 
course, more attractive to work without a special field, but it also leads to complications 
since many of the techniques to evolve inhonrogeneities are not available. Nevertheless, recent 
results EH EH suggest that this phase can generate a nearly scale invariant spectrum which 
is consistent with observations but also provides characteristic signatures compared to other 
inflation models. However, this phase alone cannot lead to a large enough universe (unless one 
assumes the parameter j to be unnaturally large) such that we need a second phase of accelerated 
expansion whose properties are not restricted so tightly since it would not give rise to observable 
anisotropies. 

Such a second phase of accelerated expansion also follows naturally from loop quantum 
cosmology with matter fields: The effective Klein-Gordon equation (HU has a </>-term which 
changes sign at small scales. This means that the usual friction in an expanding universe turns 
into antifriction very early on m ■ Matter fields are then driven away from their minima in 
the first inflationary phase and, after this phase stops and the classical equations become valid, 
slowly roll down their potentials (Fig. El- As usually, this will lead to a second (or more) 
inflationary phase which makes the universe big. 

This effect also applies if we do have an inflaton field. It will then be driven to large values 
in the loop inflationary phase, providing the necessary large initial values for the phase of slow- 
roll inflation. Moreover, since around the turning point of the scalar slow-roll conditions are 
violated, there are deviations at very early stages of the slow-roll phase which can explain the 
observed loss of power on large scales of the anisotropy spectrum EZ1EH1 










Figure 2. During the modified phase, a(t) is accelerated and cj) moves up its potential (quadratic 
in this example) if it starts in a minimum. At the peak value of the effective density, indicated 
by the dashed lines, this first phase of inflation stops, but there will be a second phase (right 
hand side) when cj) rolls down its potential before it oscillates around the minimum. Left hand 
side plots are in slow motion: each tic mark on the right t-axis stands at an increase of t by 100. 
The upper right data are rescaled so as to fit on the same plot. Units for a and </> are Planck 
units, and parameters are not necessarily realistic but chosen for plotting purposes. 


5.2. Bounces 

Both effects, the modified matter energy in the Friedmann equation and the antifriction term 
in the Klein-Gordon equation, can also lead to bounces in systems which would be singular 
classically. This provides intuitive explanations for the absence of singularities mmm and 
can be used for the construction of universe models |I21EE31EB| . 

For a bounce we need a = 0, a > 0 which is possible only if we have a negative contribution 
to the matter energy in the Friedmann equation. This can come from a curvature term or 
from a negative potential. Classically, the second condition would then be impossible to satisfy 
generically as a consequence of the singularity theorems. In the isotropic case with a scalar, this 
can also be seen from the Raychaudhuri equation m whose right hand side is negative in both 
cases: it is strictly negative with a negative potential and, if we have a curvature term, cj) will 
diverge at small a and dominate over the potential. However, when the modification becomes 
effective, the 0 2 -term in the Raychaudhuri equation can become positive m- Moreover, due to 
antifriction </> will not diverge in the closed case such that the potential can generically lead to 
a positive a m 

5.3. Quantum degrees of freedom 

The modifications used so far only relied on non-perturbative corrections coming from the 
finiteness of density operators. (In this context also perturbative corrections from the inverse 
scale factor appear above the peak, but so far they do not seem significant for cosmology 0510 
In addition, there are also perturbative corrections which are analogous to higher order or 
derivative terms in an effective action. Since the methods underlying loop quantum gravity are 
canonical, deriving an effective action is not possible in a direct manner. Nevertheless, there 



















Figure 3. Contour lines of a wave packet starting at the bottom and moving toward the classical 
singularity (dotted line) where it bounces off, compared to a solution to the effective classical 
equation (thick dashed line) derived under the assumption of a Gaussian. After the bounce, the 
wave function departs rapidly from a Gaussian and deviates from the effective classical solution. 


are methods to derive the corresponding terms in Hamiltonian equations of motion mm, and 
the appearance of quantum degrees of freedom, as would be the case for effective actions with 
higher derivative terms, can also be seen here. 

Since there are usually many different correction terms, it is not always easy to tell which one 
is dominant, if any at all. This is different from the non-perturbative effects in the density which 
can be studied in isolation by choosing a large value for the ambiguity parameter j. This is not 
always possible for other corrections, but one can test them numerically as demonstrated in m 
where a numerical evolution of the wave function under the Hamiltonian constraint in coordinate 
time has been compared with solutions to the effective classical equations with non-perturbative 
as well as perturbative corrections. One example is a term resulting from the spread of the wave 
packet which can give another explanation for bounces (Fig. OJ). 

A different method to derive effective classical equations is obtained from WKB techniques 
mm- Here, it is more difficult to deal with new degrees of freedom arising in higher order 
approximations. 

6. Less symmetric models 

All the techniques described so far for the isotropic model are now also available for homogeneous 
but not necessarily isotropic models |5T , 52] . Again, this leads to a picture at small scales very 
different from the classical one. For instance, the Bianchi IX model looses its classical chaos 
when modifications from effective densities in its curvature potential are taken into account 

jSSJEH- 














This also allows conclusions for the general situation of inhomogeneous singularities where the 
Bianchi IX model plays an important role in the BKL scenario m • Here, space is imagined as 
being composed of almost homogeneous patches each of which evolves according to the Bianchi 
IX model. Thus, if the Bianchi IX model becomes non-singular and the BKL picture remains 
valid, one can expect that singularities in general are removed. However, these assumption of a 
BKL picture available even at the quantum level is very strong such that a definite conclusion 
can only be reached by studying inhomogeneous midi-superspace models and eventually the full 
theory. Some inhomogeneities are now under control, most importantly in spherical symmetry 
EE 2L SB] which also allows conclusions for black holes 122 > SDi- The singularity issue, however, 
is more complicated in such a situation and remains to be resolved. 

Coming back to the BKL picture, we can see that not only the structure of classical 
singularities but also the approach to them is changed dramatically in effective loop cosmology. 
The classical Bianchi IX chaos implies that patches in the BKL picture have to be subdivided 
rapidly if the almost homogeneous approximation is to be maintained. This goes on without 
limit to the patch size which implies unlimited fragmentation and a complicated initial geometry 
classically. Without the chaos in the loop model, the fragmentation stops eventually giving rise 
to a minimal patch size. As can be seen [S3, the minimal patch size is not smaller than the 
scale of discreteness in loop quantum gravity thus providing a consistent picture. 

7. Conclusions 

We have reviewed the basic ingredients of loop quantum cosmology and discussed 
phenomenological consequences. Here we focused mainly on new developments after the last 
report m ■ These are in the context of structure formation from loop effects even without an 
inflaton, new terms in effective classical equations, and better techniques for inhomogeneous 
models. In particular the latter will be developed further in the near future which will not only 
bring new ingredients for cosmological investigations but also new applications in the context of 
black holes. 
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